Topological Gauge Structure and Phase Diagram of a Doped Antiferromagnet 
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We show that a topological gauge structure in an effective description of the t-J model gives rise 
to a global phase diagram of antiferromagnetic (AF) and superconducting (SC) phases in a weakly 
doped regime. Dual confinement and deconfinement of holons and spinons play essential roles here, 
with a quantum critical point at a doping concentration x c ~ 0.043. The complex experimental 
phase diagram at low doping is well described within such a framework. 
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Introduction. Cuprate superconductors have shown 
different ordering tendencies as the hole concentration 
x varies 

Q S H 3 The parent state is a Mott insula- 
tor with an AF long range order (AFLRO). Hole doping 
leads to the disappearance of AFLRO at xq ~ 0.02. In 
a lightly doped region, x < x < x c ~ 0.05, the low 
temperature (T) state is a cluster spin-glass phase with 
localized holes. At x > x c , the holes are delocalized and 
the ground state becomes a d-wave SC phase. 

In literature, based on the three-band model, dipole 
defects induced by holes have been conjectured 0, El E} 
as responsible for the destruction of the AFLRO as well 
as the spin-glass phase @ . A different cause for destroy- 
ing the AFLRO has been attributed to holes dressed 
with vortices 01- But a systematic model study is still 
lacking, which prevents us to fully understand the evolu- 
tion of the AF and SC phases at low doping. 

Theoretically, the single-band t — J model has been 
widely used to describe electronic properties in cuprates. 
But it is notably difficult to conduct a reliable inves- 
tigation continuously going from half-filling to a super- 
conducting phase. Thus, a tractable model of the doped 
AF-Mott insulator is called for. Such a model should well 
describes antiferromagnetism at half-filling, on one hand, 
and give rise to a d-wave superconductivity at large dop- 
ing, on the other hand, while the Hilbert space remains 
restricted with x as the natural charge carrier concentra- 
tion. 

An effective description which satisfies the above cri- 
teria has been derived |Tlj from the t — J Hamiltonian 
based on the bosonic resonating-valence-bond (b-RVB) 
pairing [HEH. It is given by H string = H h + H s with 

H h = -t h J2(e zA ^~ l ^)hlh J+ H.c. (1) 
H s = -J a y £(e^)blb]_ tr + H.c (2) 

Here hi and bi„ are bosonic holon and spinon opera- 
tors, respectively. At half- filling, Hh is absent while 
H s reduces to the Schwinger-boson mean-field Hamilto- 
nian 0| (Ay = 0), which correctly characterizes AF 
correlations. In particular, a spinon Bose condensation, 



< bi a >7^ 0, leads to an AFLRO. On the other hand, 
when holons are Bose condensed at large doping, the 
ground state becomes a d-wave SC [H EH- Without 
considering the spinon condensation, such a phase would 
extrapolate at x — > with T c — > 0. But at low doping, 
spins can still remain ordered. It thus may push the SC 
phase boundary to a finite x c , leaving a region for a rich 
competing phenomenon. 

In this Letter, we explore such a regime lying between 
the half-filling antiferromagnetism and the superconduct- 
ing phase. We find that the theoretical phase diagram 
bears striking similarities to the experimental one, and 
is controlled by a dual confinement-deconfinement pro- 
cedure determined by the topological gauge structure in 
Q and J2J. 

Topological gauge structure. The nontriviality of 



H 



string 



arises from the link variables, Af 



and Ay, 



the doped case. These link variables satisfy topological 



conditions: £ C A£ } . = ±7r Eie c (™ 
±tt£ 



)and £ c Ay 



lec'H 



for a closed loop c (here nf a and nf denote 
spinon and holon number operators, respectively). So Afj 
reflects the frustrations of spin background on the kinetic 
energy of charge degrees of freedom, while A\^ represents 
the influence of charge part on spin degrees of freedom. 
[tp-j in describes a uniform ir flux per plaquette with 

It has been previously shown |l4| that in the supercon- 
ducting state with a holon Bose condensation, spinons are 
"confined" , due to the gauge field Af,- , to form integer 
spin excitations and nodal quasiparticles emerge as re- 
combined holon-spinon composites. The deconfinement 
of spinon pairs occurs at T > T c , which is responsible 
for destroying the SC phase coherence. In the following, 
we show that when spinons are condensed, holons will be 
"confined" too, due to the gauge field Ay . 

Holon confinement. If spinons are condensed, it is 
straightforward to observe that a holon will cost a loga- 
rithmically divergent energy in H s . Using an expansion 
e i<M « ~ 1 + ioAy - (Ay) 2 /2, the energy cost of a holon 
in H , can be estimated as 



AE S ~ J s < &t 



>< &t > £(4) 2 

(«•> 
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FIG. 1: In the spinon Bose condensed phase, a holon-meron is 
confined with an antimeron to form a hole-dipole composite, 
which carries a charge +e and spin 1/2. The arrows denote 
n l a(-l)'{S l > 



where L is the size of the sample and a is lattice constant. 
Such an infinite energy can be physically understood as 
follows. In terms of the spin flip operator 



(4) 



the spin polarization (S^) is twisted away from the Neel 
order (— l) 1 < b\^ X bii > by an angle <&f = J^i^i Imln 
(zi — zi) nf, which has a "meron" configuration as ^ — > 
<&f ± 2ir by going around a holon once. In other words, 
a holon defined here is a topological object carrying a 
spin- meron twist. So a single holon will not appear in 
the low-energy spectrum and must be confined in the spin 
ordered phase. 

It is easy to show that two holons will repulse each 
other logarithmically in as they are merons of the 
same topological charge The only way to realize a 
finite-energy hole object in the present framework is for 
each holon-meron to "nucleate" an anti-meron from the 
vacuum. Such an anti-meron corresponds to a twist 



< b la >^< b la > exp 



i-Imln(zi - z ) 



(5) 



with zq denoting its position. Then, the meron twist $^ 
in (Sf ) or the link variable in AE S will be canceled 
out by the corresponding anti-merons if zo's approaches 
to the holon positions. Each doped hole will then behave 
like a dipole composed of a confined pair of meron (holon) 
and anti-meron as shown in Fig. 1. 

The total energy of a dipole can be generally written 
bypl 



Edipole = -E C ore + E cwc + V, (6) 

where the infinite energy in Eq. J3} is replaced by a finite 
one 



in which <5R denotes the spatial separation of two poles of 
the dipole. Here E^ ore and -E™ re represent, respectively, 
the core energies of a holon-meron and an anti-meron 
[-E^ore also includes the hopping energy from JJJ]. 

Therefore, a holon must be confined to an anti-meron 
to form a hole-dipole object in the spin ordered phase. 
Note that a dipolar result was first found by Shraiman 
and Siggia 16j based on a semi-classical treatment of the 
t — J model. But there is an important difference: in 
the present case the hole sits at a pole instead of the 
center of a dipole as shown in Fig. 1. A more detailed 
and elaborate discussion on physical properties of hole- 
dipoles in the spinon condensed phase is given elsewhere 



Quantum critical point (QCP) x c . The holon confine- 
ment holds only at low doping. With the increase of x, 
as more and more hole-dipoles are present at T — 0, the 
"confining" potential Q can get fully screened for large 
pairs of meron and anti-meron, leading to a topological 
transition at x — x c , from dipoles to free merons in a 
fashion of the Kosterlitz-Thouless (KT) transition |l7|. 

We shall employ a standard KT renormalization group 
(RG) method to calculate x c . For this purpose we rewrite 
interaction V(r) in © as 2Trf3~ 1 K\n(r/a), where the re- 
duced spin stiffness K(a) cc J s [3 (/3 = 1/T). The proba- 
bility of creating a meron-antimeron pair with two poles 
separated by a distance a is given by the pair fugac- 
ity y 2 (a). Note that in the conventional KT theory, 



V 2 {a) 



but in the present case the 



dipole number is fixed at x per site and thus the initial 
y 2 (a) must be adjusted accordingly (see below). 

In the RG scheme, small pairs of sizes within r and 
r+dr are integrated out starting from the lattice constant 
r = a. The renormalization effect is then represented by 
renormalized quantities X(r) = -j^yj and y 2 (r), which 
satisfy the famous recursion relations [TtI ITsj 



dy 2 /dl = 2(2 -^)y 2 



3„,2 



dX/dl = \-K A y 



(8) 
(9) 



where r — ae l . What makes the present approach differ- 
ent from the conventional KT theory is the presence of a 
finite density of the hole-dipoles even at T = as pointed 

out above. Here, by noting v J{' d 2 r as the areal density 
of pairs of sizes between r and r 
following constraint 



dr 01 j we have the 



x/a 2 



dr 2irr 



y 2 {r) 



2Tr 2 a 2 



(10) 



V(\SR\)=AE.~ J s m(|JR|/a), 



(7) 



In obtaining the second line the recursion relations (JSJ 
and © are used. 

The RG flow diagram of iJHJ and © is well known: 
the two basins of attraction are separated by the initial 
values which flow to X* = it/ 2 and y* = in the limit 



3 



I -► oo. For T — > 0, X(/ 
RG flows is given by 



0) — > 0, the separatrix of the £/a 



exp(|^). The Neel temperature Tn{x) is roughly 



I = 



x 



dX' 



o 4(X' - 7r/2) - 2tt ln(2X'/7r) : 



(11) 



and the critical hole density can be numerically deter- 
mined in terms of (|10f> and as 



0.84 
2^" 



0.043. 



(12) 



Therefore, a QCP is found at x c where hole-dipoles dis- 
solve into holon-merons and antimerons. At x < x c , 
since the fugacity y is always renormalized to zero, each 
holon-meron has to be bound to an immobile antimeron, 
implying that the doped holes should be self-trapped in 
space |l5j . 

Disappearance of AFLRO. Even in the confined phase, 
x < x c , the system is not necessarily always AF ordered. 
The presence of hole-dipole can lead to the destruction of 
the true AFLRO before reaching x c . The basic physics 
reason is due to the fact that dipoles have a long-range 
effect (-) on the distortion of the magnetization direc- 
tion. 

Let us introduce a unit Neel vector n; oc (— l) 1 (Si) . 
Define nf + in\ = e l< A>+ l< A° 5 with no = (cos0o,sin0o) 
as the global AF magnetization direction. Then a hole- 
dipole centered at the origin [Fig. 1] will give rise to 



(13) 



at |rj| >> |p|, with p = — zx^R . 

The twist of the Neer vector due to a dipole at the orig- 
inal is given by <5n»= n;— no ~ va.<j)i, with the unit vector 
m = (— sin 0o, cos^o)- In the continuum limit, V 2 n(r) = 
27rm (p ■ V) <5(r). The multi-dipole solution can be gen- 
erally written as Vn M (r) = 27rX); Pl${ r — r l) m i + S± 
where I denotes the index of hole-dipoles and V • g^ = 0. 
The transverse component g^ will have no effect in the 
nonlinear sigma model 7] and we may only focus on 
the longitudinal part of Vn M (r) below. At low tem- 
perature, we may assume that the hole-dipoles are lo- 
calized and treat all the variables, ni;, p;, and rj, as 
quenched. Defining the quenched average (• • •) , and 



using (mfmY, 



77= (\SK\ 2 ) q /a 2 , we get 



l/28^ v 5 u >, (pIp 3 ^ = 5 lj 5 ll >r]a 2 /2 with 



(*»'<^nV)> g = W(r-r') (U, 

in which v = Ax, with A — Tr 2 r). 

The RG study of the non-linear sigma model with 
quenched random dipole moments has been given 
within a one-loop approximation. Even though the ori- 
gin of the dipole moments is different, once l|14|l is de- 
termined, these results can be directly applied here. The 
AF correlation length £ has been obtained at low T as 



given by the solution of a£ 2 ~ a 2 , where a ~ 10 -5 Q|, 
representing the effect of the interlayer coupling J±/J. 
Then the critical doping xq at which the AFLRO disap- 
pears can be estimated by xq = — %4j- - ■ In order to get 
the experimental value xq ~ 0.02, it has been assumed 
A ~ 20 in Ref. 0. In the present self-consistent 
calculation leads to rj = 1.23 and xq ~ 0.03, determined 
by the KT theory based on © and ©. On the other 
hand, at x — > where v <C t = Tm / p s (p s ~ 0.176 J is the 
spin stiffness [l3|), f/a ~ exp(§J[l-(l-f ) 3 ]) 7] and one 
obtains T N (x) « T N (0) - Axp s , with T N (0) = - fep,. 
The plot of Tjy as a function of x is shown in Fig. 2. 

In Fig. 2, a characteristic temperature Tf in the 
AFLRO phase is also shown by the dotted curve, which 
represents the fact that although the holes are all local- 
ized, the directions of their dipole moments can still ro- 
tate freely to reach annealed equilibrium above Tf. The 
dipole-dipole interaction causes an energy difference of 
two dipoles, from parallel to perpendicular in their rela- 
tive moment alignment, is proportional to 1/r 2 (r is the 
spatial separation between them). Associating r with the 
average hole-hole distance: r = a/y/x, the interaction en- 
erg y th en scales linearly with x, such that Tf ~ 1/r 2 ~ x 
[H . In the region xq < x < x c , the AFLRO is 
destroyed and the AF orders are limited mainly by fi- 
nite size effects, where the size of the AF domains is 
determined by hole concentration, £ ~ a/ ' yfx. The spin- 
glass freezing temperature is then expected to vary as 
^9 ~ £ 2 ~ I/ 2 -- Such a phase has been known as a cluster 



spin glass |8| • Below the temperature T g the hole-dipolar 
configurations form a glass and their dynamics strongly 
slows down. 

Superconducting phase. At x > x c , holons are decon- 
fined and free, and thus will experience Bose condensa- 
tion at low T, giving rise |lll Il4| to a d-wave supercon- 
ducting ground state. 

We have shown that the holons are confined to form 
immobile dipoles at x < x c , instead of being Bose con- 
densed. At x — x c , the concentration of anti-merons is 
x c , while at x > x c „ for each additional holon, there 
will be no more anti-meron to be created so that the 
concentration of anti-merons is roughly fixed at x c until 
<| b ia |>= 0. Thus, in an overlap regime of <| b ia |>^ 
and < hi >^ at x > x c , A^ in H s should be re- 
placed by a A!£j due to the presence of antimerons, and 

its strength B h = irx e ff/a 2 is now controlled by an ef- 
fective concentration instead of x itself. 
Eventually at higher doping, <| 6; CT |>= with no more 
(anti)merons, one will recover x from x e ff as the param- 
eter representing the doping effect on the SC properties 



At small doping, the continuum version of ^ is a 
CP 1 [2(j model, which leads to the Klein-Gordon equa- 

, 2 

tion: 

ha) 

tices 



(a, 



iA' l a 



- [E/c s ) 2 z a , where z a = \(b Aa + 

b Aa and bg a are spinons on A and on B sublat- 
and c s is the effective spin-wave velocity (c s ~ 



4 



T/J 




FIG. 2: Phase diagram at low doping x : a dual confinement- 
deconfmement occurs at a quantum critical point x c ~ 0.043. 
The Neel temperature Tjv vanishes at xq ~ 0.03. T/ and T 9 
denote characteristic spin freezing temperatures, and T c is the 
superconducting transition temperature. E g is a characteris- 
tic spin energy (see text). 



This result is in contrast to E s oc xJ (x — > 0) obtained at 
<\bicr |>= 0| ■ The characteristic spin energy is then 
given by E g ~ 2£J S . It has been previously established 
that the superconducting transition occurs when 
spinons become deconfined at a finite temperature, which 
is determined at T c ~ with c ~ 4, as shown in Fig. 2. 

In conclusion, the low-temperature phase diagram for 
a doped Mott antiferromagnet described by (JIJ and © 
is basically determined by its intrinsic topological gauge 
structure. A quantum critical point at x c ~ 0.043 is 
found, below and above which, dual confinement and de- 
confinement take place at T = 0, leading to a systematic 
evolution from the antiferromagnetic to superconduct- 
ing phases as a function of x. The complex experimen- 
tal phase diagram in the weakly doped cuprates may be 
understood within such a framework. For example, the 
dipolar effect of doped holes is responsible for a vanish- 
ing T/v at xq ~ 0.03 and a cluster spin-glass phase at 
xq < x < x c . The superconducting state sets in at x > x c , 
whose phase coherence is destroyed by a deconfinement 
of spinons at T > T c [23l l24j . Due to the space limit, 
we have not discussed the possible stripe instability in 
the same model, which is explored at low doping as a 
competing phenomenon elsewhere |Tflj . 



IMaJ [13). The energy spectrum is [13, E n>m = 

±c s \J 1 (n + 1/2+ | m | ^fm)B h and the wave functions are 

<,m ~ P |m| e iro< ^F(-ri,| m | +1, a 2 p 2 ) exp(-aV/2), 
where F is the hypergeometric function, m is the eigen- 
value of the angular momentum, and n is the eigen- value 
of the harmonic oscillation level number, a 2 = B h /2, 
p = \J x 1 + y 2 . The energy gap between first excited 
state and the ground state is 

E s = E 1 -E a ^ UaJy/x^j. (15) 
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